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Abstract 


We  used  the  GISS/Harvard/UCI  3-D  CTM  on  the  AER  computer  platform  to 
perform  a  number  of  simulations  to  obtain  the  surface  area  from  accumulation  of  AI2O3 
particulates  from  solid  rocket  motors  (SRM),  orbital  debris,  and  meteorites.  In  the 
calculation,  the  initial  size  distribution  of  the  particulates  emitted  by  SRM  is  represented 
by  a  tri-modal  distribution  with  bulk  density  of  1.7  gm/cm  .  It  is  assumed  that  particles 
do  not  interact  with  each  other  so  that  particles  in  each  size  bin  will  evolve  independently. 
Particles  from  orbital  debris  and  meteorites  are  assumed  to  be  1  pm  size  particles  with 
bulk  density  of  4.7  gm/cm  and  2.0  gm/cm  ,  respectively.  Apart  from  the  large-scale 
transport,  particle  distributions  are  also  affected  by  sedimentation.  In  addition,  we 
consider  the  effect  of  removal  of  the  particles  by  collision  with  the  background  sulfate 
layer. 


The  calculated  surface  area  of  AI2O3  particulates  from  SRM,  with  an  input  of 
1120  tons/yr,  is  between  1-4x10"  pnr/cm  in  the  lower  stratosphere  with  the  largest 
values  occurring  at  northern  high  latitudes.  The  assumed  input  from  orbital  debris  and 
meteorites  are  much  smaller  at  10  tons/yr  and  30  ton/yrs  respectively.  The  calculated 
surface  areas  are  also  much  smaller  at  10  pm "/cm  and  10"  pnr/cm  respectively. 

2  3 

These  values  are  to  be  compared  with  the  sulfate  surface  area  of  0.5- 1.0  pnr/cm  in  the 
lower  stratosphere  during  non-volcanic  conditions,  and  value  of  20  pnr/cm  one  year 
after  the  Pinatubo  eruption. 

The  accumulation  of  AI2O3  particulates  in  the  atmosphere  may  affect  ozone  via 
the  heterogeneous  reaction  C10N02  +HC1  — >  HN03  +  Cl2  that  converts  CIONO2  and 

HC1,  chlorine  reservoir  species,  into  the  more  active  form  that  will  deplete  ozone  in  the 
presence  of  sunlight.  However,  due  to  the  very  small  A1203  surface  areas  calculated, 
ozone  depletion  on  the  global  scale  is  very  small.  The  impact  on  the  stratospheric  sulfate 
aerosol  layer  is  likely  to  be  small  over  most  of  the  stratosphere,  but  the  impact  cannot  be 
evaluated  with  accuracy  at  this  time. 
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Section  1:  Introduction 


Description  of  the  global  dispersions  of  materials  emitted  by  rockets  involves 
simulations  at  different  spatial  and  temporal  scales.  The  primary  interest  in  this  study  is 
the  global  scale.  For  the  purpose  of  this  study,  we  will  assume  that  emitted  materials 
from  the  rocket  are  dispersed  instantaneously  and  uniformly  through  the  model  grid  box 
(see  below  for  the  size  of  the  model  grid-box). 

Calculations  in  this  report  were  performed  using  the  21  layer  GISS/Harvard/UCI 
Three-Dimensional  Chemistry  Transport  Model  (3-D  CTM).  This  model  is  upgraded 
from  the  9  layer  Harvard/GISS  3-D  CTM  for  the  troposphere,  which  has  been 
documented  by  Prather  et  al.,  [ 1987].  The  21  layer  model  has  been  described  in  detail  in 
Prather  et  al.,  [ 1990a].  Both  models  have  been  used  to  study  a  variety  of  problems  in  the 
atmosphere,  which  include  CFC  simulations  (Prather  et  al.,  [1987]),  the  effects  of  debris 
from  meteors  on  the  Antarctic  ozone  (Prather  and  Rodriguez  [1988]),  the  dynamical 
dilution  of  the  Antarctic  ozone  hole  (Prather  et  al.,  [1990a]),  the  Space  Shuttle's  impact 
on  the  stratosphere  (Prather  et  al.,  [ 1990b]),  the  trend  and  annual  cycle  in  stratospheric 
CO2  (Hall  and  Prather  [1993]),  the  seasonal  evolution  of  N2O,  O3,  and  CO2  (Hall  and 
Prather,  1995)  and  trace-tracer  correlation  in  the  stratosphere  (Avallone  and  Prather 
[1997]). 


1.1  Model  Resolutions  and  Winds 
A.  Model  resolution  and  range 

This  is  a  global  model  with  a  resolution  of  7.83°  (latitude)  x  10°  (longitude)  x  21 
layers  (altitude).  In  the  longitudinal  direction,  there  are  36  boxes  (1=1  to  36)  with  10 
degree  longitudinal  width.  The  1=1  box  is  centered  on  the  date  line  (180E  =  180W)  and  I 
increases  eastward.  In  the  latitudinal  direction,  there  are  24  boxes  (J=l  to  24)  from  pole  to 
pole  with  the  J=1  box  at  the  South  Pole.  The  7.83°  width  is  equal  to  180723  to  get  an 
'over-the-pole'  box  of  the  same  length  as  the  other  boxes.  In  the  vertical  direction,  there 
are  21  layers  (L=l  to  21).  The  bottom  9-layers  (L=l  is  near  the  surface)  are  pure  <7  layers 
with  the  pressure  at  the  edge  defined  by: 


-3- 


p(edge)  =  <7  *  (Psurf  -  100  mbar)  +  100  mbar 


with  <7  =  1.00,  0.973,  0.938,  0.887,  0.788,  0.638,  0.460,  0.278,  0.117,  0.0.  The 
top  12-layers  (from  100  mbar  up)  are  at  fixed  pressure  levels,  with 


p(edge)  =  100.0,  46.41,  21.54,  9.998,  4.641,  2.155,  1.000,  0.464,  0.216, 
0.100,  0.0466,  0.0215,  0.0022  mbar. 


B.  Transport  numerics 

The  Prather  advection  scheme  (Prather  [1986]),  which  conserves  second 
moments  of  tracer  distribution,  is  used  in  the  model.  Although  this  scheme  requires  more 
CPU  time  and  more  storage  space  to  run,  previous  work  (see  Shia  et  al,  [ 1990] )  shows  it 
does  a  good  job  for  simulating  advection  of  tracers  with  large  gradients.  The  only 
(explicit)  diffusion  included  in  the  model  is  the  horizontal  diffusion  caused  by 
convection. 

C.  Meteorological  fields 

The  3-D  CTM  is  driven  by  the  meteorological  fields  from  the  GISS-II  general 
circulation  model  (GCM)  with  23-layers,  GCMAM,  that  includes  basic  GISS-II  GCM 
physics  for  the  9-layer  model  with  gravity  wave  drag  for  the  stratosphere  and  other 
improvements  (Rind  et  al.,  [ 1990],  Prather  et  al.,  [ 1993 ] ).  The  top  3  layers  have  been 
combined  into  a  single  layer  to  drive  the  21-layer  CTM.  We  have  received  the 
meteorological  fields  generated  by  GCMAM  every  8  hours  for  365  days.  The  data  are 
used  repeatedly  every  year. 


1.2  Performance  on  the  AER  platform 

We  received  from  Prof.  Prather  of  UCI  a  CD-ROM,  which  includes  the 
FORTRAN  source  code  for  the  21  layer  GISS/Harvard/UCI  3-D  CTM  and  also  the  inputs 
and  outputs  for  two  sample  runs  for  testing  (Prather,  private  communication,  [1997]). 
The  code  has  been  ported  on  to  a  SGI  Power  Challenge  machine  at  AER.  First,  the  two 
sample  runs  were  repeated.  Both  runs  simulate  N2O  and  NOy  distributions  in  the 
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atmosphere,  one  for  7  years  and  one  for  1  week.  The  shorter  run  has  more  detailed 
diagnostic  output,  including  the  global  budget,  the  tendency,  the  concentrations  of  N2O 
and  NOy.  Our  model  outputs  are  the  same  as  those  calculated  at  UCI,  which  came  with 
the  CD-ROM.  The  CPU  time  (one  processor)  for  the  7  year  run  on  the  SGI  machine  is 
about  9  hours. 

In  order  to  use  this  model  to  simulate  the  transport  of  an  inert  tracer,  like  the 
emissions  from  rockets  or  deorbiting  space  debris,  all  the  sources/sinks  in  the  model  have 
been  deactivated.  For  an  inert  tracer,  if  there  is  no  deposition  to  the  surface  the  total  mass 
of  the  tracer  should  be  conserved.  We  tested  the  mass  conservation  and  after  one  year 
model  run  the  total  mass  of  the  inertial  trace  is  conserved  closed  to  the  machine  accuracy 
(less  than  10' 10). 


1.3  Numerical  Experiments  Simulating  Accumulation  of  Chlorine  from  Rocket  Launches 

We  performed  two  experiments  to  simulate  the  distributions  of  chlorine  emitted 
by  rocket  launches.  The  source  term  in  our  experiment  corresponds  to  the  chlorine 
emissions  from  nine  launches  of  the  Space  Shuttle  yearly.  The  source  is  located  above 
Cape  Canaveral,  Florida  (29°N,  80°W)  and  its  vertical  distribution  is  taken  from  Prather 
et  al.,  [1990b].  The  inputs  are  concentrated  in  northern  mid-latitudes  and  the 
magnitudes,  in  mixing  ratio  unit,  are  largest  around  30-40  km  region.  In  the  simulation, 
the  mixing  ratio  of  the  tracers  in  the  lowest  3  layers  is  continuously  set  to  zero  to  simulate 
efficient  washout  to  the  ground. 

In  the  first  experiment,  the  emissions  are  put  into  the  model  every  40  days  starting 
at  day  40  every  year  to  represent  each  individual  launch.  The  input  of  Cl  into  the 
stratosphere  is  68  ton  for  every  launch  and  the  total  inputs  of  Cl  every  year  is  68  ton  x 
9  times/year  =  612  ton/year.  After  10  years,  the  tracer  distribution  reaches  an  annual 
repeating  steady  state.  The  total  calculated  Cl  content  in  the  atmosphere  is  1413  ton, 
which  means  that  the  residence  time  of  Cl  from  Space  Shuttle  launches  is  about 
1413[ton]  /  612[ton/year]  =  2.3  years.  The  zonally  averaged  Cl  mixing  ratio  distribution 
in  the  steady  state  for  4  seasons  is  plotted  in  Figure  1A.  Figure  1A  shows  the  dispersion 
of  the  input  Cl  downward  and  upward  from  this  region  and  also  across  the  equator  into 
the  Southern  Hemisphere.  To  assess  the  effect  of  the  intermittent  source,  we  performed  a 
second  experiment  in  which  the  emissions  from  the  same  nine  launches  were  spread 
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uniformly  through  the  year.  The  results  from  this  experiment  are  shown  in  Figure  IB. 
Comparison  of  Figures  1A  and  IB  indicates  that  there  is  no  significant  difference 
between  the  calculated  distributions.  Since  the  exact  timing  of  the  launches  are  not  fixed, 
we  will  use  the  continuous  source  function  in  the  rest  of  the  calculations 

The  seasonal  variability  of  the  transport  is  demonstrated  in  Figures  1A  and  IB. 
Due  to  the  weaker  transport  in  later  summer  and  early  fall  in  the  Northern  Hemisphere, 
the  input  Cl  tends  to  accumulate  the  most  in  summer  and  fall  in  the  source  region,  a 
maximum  of  the  Cl  mixing  ratio  of  1 1  pptv.  Results  from  Prather  el  al.,  [ 1990b]  is 
shown  in  figure  1C  for  comparison.  The  differences  between  Figures  IB  and  1C  are 
caused  by  the  differences  in  magnitudes  and  locations  of  source  input  (nine  shuttle 
launches  versus  nine  shuttle  launches  plus  six  Titan  launches). 
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Section  2:  Overall  Approach 


In  this  project,  we  will  treat  the  particulate  from  the  various  sources  as  non¬ 
interacting  particles.  It  is  important  to  list  the  assumptions,  along  with  the  rationales,  that 
will  be  made  in  this  study. 

The  main  purpose  of  the  study  is  to  provide  an  estimate  of  the  effect  on 
stratospheric  ozone  from  reactions  that  may  occur  on  the  surfaces  of  the  particles  from 
rocket  exhaust  and  space  debris.  Since  the  reaction  rates  are  proportional  to  the  surface 
area,  we  can  obtain  an  upper  limit  on  the  effects  by  maximizing  the  surface  area.  In  the 
atmosphere,  particles  may  coagulate  to  form  larger  particles.  Coagulation  will  reduce  the 
surface  area  (larger  particles  have  a  smaller  surface  to  volume  ratio).  Larger  particles 
have  larger  sedimentation  velocity  and  are  removed  more  rapidly  from  the  atmosphere. 
Thus  ignoring  coagulation  in  the  simulation  will  maximize  both  the  surface  area  density 
and  mass  density  of  the  particles  in  the  atmosphere.  Although  there  are  indications  that 
H2SO4  may  not  readily  condense  on  ALO3  particles,  we  feel  that  it  is  nonetheless  useful 
to  estimate  the  effect  of  the  particles  as  condensation  nuclei  for  sulfate  aerosol.  Again,  in 
this  case,  smaller  particles  are  more  efficient  as  condensation  nuclei.  For  these  reasons, 
we  decided  that  we  will  treat  the  emitted  particles  as  non-interacting  and  ignore 
coagulation  in  the  calculations.  Once  the  emitted  particles  are  dispersed  to  the  size  of  the 
model  grid  box,  number  densities  will  be  low  enough  that  coagulation  between  A1 
particles  will  not  be  an  important  process  anyway. 


2.1  Sedimentation  Velocity 


Sedimentation  is  an  important  process  for  atmospheric  particles,  affecting  their 
residence  time  and  vertical  distribution.  Our  treatment  of  gravitational  sedimentation 
follows  Kasten  [ 1968].  The  falling  speed  (wse(j)  of  spherical  particles  due  to  gravity  is  a 
function  of  particle  density  ( phulk),  radius  (D),  air  viscosity  ( rj)  and  mean  free  path  (A): 


|  Phu,kD2g  | 

r,+4f 

9 

L  D\ 

1.249  +0.42  exp 


0.87D 

A 


(1) 


Air  viscosity  is  a  function  of  ambient  temperature,  while  mean  free  path  is  a  function  of 
both  temperature  and  atmospheric  pressure.  Particles  between  0.02  pm  and  0.1  pm  radius 
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do  not  settle  appreciable  in  the  lower  stratosphere,  but  are  influenced  by  gravitational 
sedimentation  at  higher  altitudes.  Particles  greater  than  0.1  pm  radius  are  influenced  by 
sedimentation  at  all  altitudes,  and  particles  greater  than  1.0  pm  radius  are  rapidly 
removed  from  the  stratosphere.  The  sedimentation  velocity,  defined  in  Equation  (1),  is 
plotted  as  a  function  of  the  particle  size  in  Figure  2.  The  three  curves  represent  the 
sedimentation  velocity  at  20°  N  around  20  km,  30  km,  and  40  km  respectively.  At  the 
same  altitude,  bigger  particles  have  larger  sedimentation  velocity  because,  relatively 
speaking,  the  air  resistance  is  smaller  for  the  bigger  particles.  For  particles  with  the  same 
size,  the  sedimentation  velocity  is  larger  at  higher  altitude  because  the  thinner  air  at 
higher  altitude  provides  less  friction  for  the  particles.  The  sedimentation  velocity  is 
almost  20  times  larger  at  40  km  than  that  at  20  km.  Typical  magnitudes  of  the  vertical 
advection  velocity  from  the  large  scale  circulation  are  0.05  km/day  at  20  km,  and  0. 1 
km/day  at  40  km.  Thus  sedimentation  is  important  only  for  particles  larger  than  0.1  pm. 


2.2  Removal  by  collision  with  sulfate  aerosol 


A  first  order  loss  rate  for  each  size  bin  of  aluminum  particles  can  be  generated 
using  the  following  formula: 


dt 


KjkNk 

k>j 


sulfate 


(2) 


where  rij  represents  the  number  density  of  A1203  particles  in  the  jth  size  bin,  N™'-fate 
represent  the  number  density  of  background  sulfate  particles  in  the  kth  size  bin,  and  Kjk 
represents  the  collision  rate  (Yue  &  Deepak  [ 1979])  between  bins  j  and  k.  The  size 
distribution  of  sulfate  particles  for  background  (non-volcanic)  conditions  is  taken  from  a 
2-D  model  calculation  by  Weisenstein  et  al.,  [1997].  We  assume  in  equation  (2)  that 
A1203  particles  only  coagulate  with  sulfate  particles  of  same  or  larger  size,  ensuring  that 
the  AFO3  particle  is  completely  coated  with  sulfate  before  it  is  removed.  By  adopting 
this  assumption,  we  excluded  the  possibility  that  the  A1203  particles  may  eventually  be 
coated  by  a  series  of  collisions  with  smaller  sulfate  particles.  Thus,  the  results  should 
represent  an  upper  limit. 


The  collision  rate  tends  to  be  largest  for  collisions  between  small  and  large 
particles,  because  small  particles  have  a  high  thermal  velocity  and  large  particles  have  a 
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large  cross-section.  Therefore,  only  the  small  AI2O3  particles  are  removed  efficiently  under 
these  assumptions.  This  characteristic  is  shown  in  Figure  3,  in  which  the  collision  removal 
rate  is  plotted  as  function  of  particle  size  at  20  km  for  the  equator,  35°  N  and  70°  N.  The 
latitude-altitude  cross-section  of  the  reciprocal  of  the  removal  rates  for  4  seasons  are  plotted 
for  particles  with  0.03  pm  radius  and  particles  with  0.13  pm  radius  in  Figures  4  and  5, 
respectively.  The  collision  removal  is  much  slower  for  particles  with  0.13  pm  radius  than 
particles  with  0.03  pm  radius  (20  to  30  times  slower).  The  time  constant  for  the  0.13  pm 
case  is  much  larger  than  the  stratosphere  residence  time  of  800  days.  Thus,  collision 
removal  should  have  little  effects  for  0.13  pm  particles. 


2.3  Model  Simulations 

We  used  the  GISS/UCI  3-D  CTM  to  simulate  the  quasi-steady  state  (annually 
repeating)  distribution  of  AI2O3  particulate  for  the  following  three  cases, 

•  Baseline:  particulate  source  transported  by  large  scale  circulation  ignoring  sedimentation, 

•  Case  A:  particulate  source  with  sedimentation, 

•  Case  B:  particulate  source  with  sedimentation  and  collision  removal  with  the 
background  sulfate  particles. 

Note  that  since  we  assume  the  A1203  particles  are  non-interacting,  the  baseline  case 
requires  only  one  model  run  for  all  sizes  of  particles.  As  there  are  no  interactions 
between  AI2O3  particles,  the  simulations  in  Case  A  and  Case  B  are  run  independently  for 
each  size  bin  and  the  results  combined  to  obtain  the  surface  area. 
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3.  Particulates  from  Solid  Rocket  Motors 


3.1.  Description  of  Source  Function 

Malcolm  Ko  (AER)  worked  with  Peter  Lohn  (TRW)  to  finalize  the  description  of 
the  source  for  solid  fuel  rockets.  Beiting  [1995,  1997 ]  reviewed  the  measurements 
through  rocket  plumes  and  determined  that  the  particles  can  be  described  by  a  tri-modal 
distribution  specified  as  follows: 

3 

nD(D)  =  nf-y>D  (3) 

1=1 

where  D  is  the  radius  of  the  particles,  n;  is  the  number  density  (particle  pm"1  cm  3)  with 
n,/n2  =  4.65  xlO4,  n2/n3  =  10.3, 


Yi  =  63.3  pm"1 

for 

y2=  3.13  pm"1 

for 

y3  =  0.80  pm  1 

for 

0.025  pm  <  D  <  0.25  pm 
0.25  pm  <  D  <  1.0  pm 
1.0  pm  <  D  <  10  pm 


If  we  define  qmass(x,t)  (gm/cm)  as  the  mass  of  aluminum  oxide  deposited  per  cm  vehicle 
track  at  location  x  and  time  t,  the  number  of  particles  introduced  into  the  model  grid  box 
per  unit  vehicle  track  q(D,x,t )  (in  units  of  particles  pm  1  cm  V 


q{D,x,t) 


(LuJAA)  3 

Pbulki*’  0  m 


(4) 


where  phuik(  x,t)  is  the  bulk  density  of  the  particles  in  gm/cc,  and  the  normalization  factor 
?7(  has  units  of  [particles  pm  cm  ].  The  nts  are  obtained  by  solving  the  following 

equations 

1  =  3  ri^NVX 

1=1 

n\  =  4.65  x104772 
n'2=  10.3 
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where  (iW),-  is  the  normalized  volume  given  by 


a  —  D=0.25/Mn 

(NV)i  =  —  e~r'DD3dD  =  1.45xl0'18  cm3-(im 

3  D=0. 025 jum 

Att  D=l.Ojum 

(NV)2  =  —  e~7lDD3clD  =  9.78xl0"14  cm3-pm 

3  D=0.25pn 

A—D=\0.0/J.m 

(NV)3  =  —  e~YlDD3dD  =  5.82xl0"n  cm3-(im 

3  D=l.0jUm 

The  solutions  are  n[=  7.85  xlO15  [particles  pm"1  cm‘3],n2=  1.69  xlO11  [particles  pm"1  cm3]  , 
and«3=  1.64  xlO10  [particles  pm"1  cm"3]. 


In  the  Beiting  [1995]  study,  it  is  assumed  that  the  particles  are  non-interacting  and 
this  tri-modal  distribution  is  preserved  as  the  plume  expands.  In  the  calculation,  the 
limited  particles  are  mixed  into  the  grid  box.  Information  from  the  TRW  team  indicates 
that  qmass(x,t)  can  be  specified  according  to  Table  1  for  Titan  IV  launches,  with  a  total 
emission  of  55  tons  per  launch.  The  same  spatial  distribution  is  used  for  shuttle  launches, 
but  with  100  tons  emitted  per  launch. 


Table  1:  Vertical  Distribution  of  the  A1203  Particulate  Source  from  each  Titan  IV  Launch 


Altitude  range 

(km) 

Mass  of  A1203  deposited 

(Ton) 

9  mass 

(gm/em) 

15-20 

11.6 

2.32xl04 

20-25 

9.2 

1.84xl04 

25-30 

9.1 

1.82xl04 

30-35 

7.4 

1.48xl04 

35-40 

7.2 

1.44xl04 

40-45 

5.7 

1.14xl04 

45-50 

4.8 

0.96xl04 
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3.2  Model  Simulations 


The  source  of  A1203  particulate  from  SRM  includes  nine  space  shuttle  launches 
every  year,  which  deposit  900  ton  of  AI2O3  particulate  into  the  stratosphere,  and  four 
Titan  launches  every  year,  which  deposit  220  ton  of  A1203  particulate  into  the 
stratosphere.  We  put  the  total  source  of  1 120  ton/yr  into  the  model  uniformly  every  day  at 
(20°N,  80°W)  with  the  vertical  distribution  listed  in  Table  1  of  Section  2.  For  this  study, 
we  decided  to  input  the  source  in  to  the  model  every  day  uniformly  is  more  suitable  than 
to  input  a  large  amount  of  source  at  a  short  period  (e.g.,  in  several  hours).  If  we  cannot 
keep  the  launch  schedule  the  same  every  year,  we  think  the  results  of  the  first  scenario  are 
closer  to  reality. 

In  this  study,  we  assume  that  the  particulate  from  the  SRM  do  not  interact  with 
each  other.  They  are  removed  by  sedimentation  into  the  troposphere  followed  by 
washout;  and/or  by  collision  with  the  background  stratospheric  sulfate  aerosols. 


Since  both  removal  rates  are  size  dependent,  we  must  perform  separate 
calculations  for  each  size.  In  the  3-D  model  simulations,  the  particles  of  A1203  from 
SRM  are  divided  into  13  bins  covering  particle  radius  from  0.025  //m  to  10  //  m .  The 
boundaries  of  bins  are  Dj+1/2  with  j  =  0-13.  Their  values  are  listed  in  the  first  column  of 
Table  2.  The  size  distribution  of  the  A1203  particulate  emitted  by  SRMs  is  described  by 
the  tri-modal  distribution  (Equation  3).  In  the  13  bins,  which  are  used  to  cover  the  size 
distribution  of  A1203  particulate  for  3-D  model  simulation,  the  first  5  bins  are  in  the  first 
mode,  bins  6-8  are  in  the  second  mode,  and  the  largest  5  bins  in  the  third  mode. 


The  total  source  (1120  ton/yr)  of  the  A1203  particles  from  SRM  can  be  expressed 
as  the  sum  of  sources  in  each  bin, 


13 

S total  ~  Sj 

j= 1 


(5) 


n  e  y‘dD3cID 

AD  j  ' 


(6) 
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where  ^  is  the  bulk  density  of  the  particles  (assumed  to  be  1.7  gm/cc),  and 


A Dj  =  (pj+l/2  -  Dj_u  2 )  is  the  range  of  the  j  th  bin, 

with  i  =1,  =4. 79x10511m'1,  yi=63. 3  pm"1,  for  7=1-5, 
with  i  =2,  «2  =  10.3pm'1,  y2=3.13  pm'1,  for  7=6-8, 
with  i  =3,  n ^  =  1.0pm'1,  Y3=0.8  pm"1,  for  7=9-13 


The  constant  C  is  obtained  by  requiring  that  Stotai  =1120  ton/yr  and  it  is  found  to  be  1.1 
x  1019  yr"1.  The  values  of  Sj  are  given  in  Table  2.  They  are  distributed  in  altitude  as 
listed  in  Table  1. 


The  model  computes  the  mass  of  the  A1203  particles  in  each  of  the  13  bin  in  each 
grid  box  at  each  time  step.  As  explained  in  Section  2.3,  we  considered  three  cases:  the 
baseline  case  with  no  sedimentation  and  no  collision  removal,  Case  A  with 
sedimentation,  and  Case  B  with  sedimentation  and  collision  removal.  Thus,  a  total  of  27 
(2  x  13  +  1)  simulations  are  performed. 

Since  we  assume  for  the  baseline  case  that  particles  of  all  sizes  are  transported  at 
the  same  rate,  the  tri-model  distribution  is  preserved  and  we  can  use  that  information  to 
resolve  the  size  within  each  bin  at  each  subsequent  time.  This  is  no  longer  possible  for 
Case  A  and  Case  B  because  the  initial  size  distribution  within  each  bin  will  be  modified. 
For  the  calculation,  we  assume  that  the  particles  within  each  bin  have  the  same  radius 
given  by 

Dj  ~(Pj- 1/2  xDj+U2^  7 


the  sedimentation  velocity  (e.g.,  1)  and  the  collision  rate  (e.g.,  2)  are  calculated  using  Dj. 


3.3  Residence  Times  of  A1203  Particulate  from  SRM 

The  residence  time  is  defined  as  the  ratio  of  the  total  burden  of  the  A1203  particles 
in  steady  state  to  the  rate  of  source  input.  Without  sedimentation  or  collision  removal, 
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the  residence  time  of  the  AI2O3  particles  is  2.7  years.  With  sedimentation,  the  lifetime  of 
large  size  particles  is  reduced  markedly.  In  contrast,  the  collision  removal  only  affects 
particles  with  small  radius.  With  these  removal  processes  operating,  the  atmospheric 
lifetimes  of  particles  are  functions  of  particle  radius.  These  are  plotted  in  Figure  6.  In  the 
baseline  simulation,  particles  of  all  sizes  have  the  same  source  and  sink.  The  only  sink  is 
removal  at  the  surface.  Therefore,  all  the  particles  have  the  same  lifetime  of  2.7  years, 
indicated  by  the  solid  line  in  Figure  6.  Note  that  this  is  longer  than  the  2.3  years 
computed  for  chlorine  (Section  1.3).  The  reason  for  this  difference  is  that  the  source 
function  for  chlorine  from  Prather  et  al.,  [ 1990b]  has  the  chlorine  emitted  at  lower 
altitudes  than  the  source  function  in  Table  1  for  AI2O3.  When  the  sedimentation  is  added 
(Case  A)  the  removal  of  large  particles  at  the  surface  is  accelerated  and  their  lifetime  is 
reduced.  The  collision  removal  only  affects  particles  of  small  size.  That  is  why  Case  A 
and  Case  B  differ  only  for  sizes  less  than  0.13  pm.  The  mass-weighted  lifetime  is  about 
0.3  years  for  both  Cases  A  and  B  since  emissions  for  particles  smaller  than  0.13  pm  are 
small. 


3.4  Surface  Area  of  AFO3  Particulate  from  SRM 


The  3-D  model  calculations  generate, mAx) ,  the  mass  of  the  A1203  particles  in 

the  j  th  bin  in  the  grid  box  (x).  The  surface  area  for  the  jth  bin  (a.(x))  can  be  calculated 
from  the  mass  and  the  number  of  AI2O3  particles  n;(x)  in  the  jth  bin. 


nj(x)  =  mj(x)/(pbulk 


a.j(x )  =  rij(x)  x  4n D]  = 


3  ntj  (x) 

Pbuik  x  Dj 


(8a) 


(8b) 


Adding  all  13  bins  together,  we  obtain  the  total  local  surface  area  density, 

13 

a(x)=  cij(x) 

1  =  1 


(9) 


The  total  surface  area  density  around  20  km  in  January  for  Case  A  is  plotted  in  Figure  7, 
which  shows  that  most  AFO3  particles  remain  in  the  Northern  Hemisphere.  The  surface 
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area  density  is  more  than  four  time  higher  in  the  northern  mid-latitudes  than  in  the 
Southern  Hemisphere.  In  the  longitudinal  direction,  the  distributions  of  the  particle  are 
quite  evenly  due  to  the  strong  zonal  wind.  The  zonally  averaged  total  surface  area 
densities  are  plotted  for  four  seasons  for  the  Baseline  Case,  Cases  A  and  Case  B  in 
Figures  8,  9  and  10,  respectively.  The  Baseline  Case  (Figure  8)  has  the  largest  surface 
area  density  and  Case  B  (Figure  10)  has  smallest  one.  This  is  because  that  the 
sedimentation  and  the  collision  removal  reduce  the  particle  density  in  the  air.  In  all  three 
cases  the  total  surface  area  density  is  smaller  in  the  tropics  and  subtropics  than  in  the  mid¬ 
latitudes  of  the  lower  stratosphere.  The  behavior  in  the  tropics  matches  the  shape  of  the 
tropopause.  The  strong  convective  mixing  keeps  the  particle  number  density  low  in  the 
upper  tropopause.  The  order  of  magnitude  of  the  total  surface  area  density  is  10"4 
( ji  m^/cm3),  which  is  about  3  orders  smaller  than  that  of  the  background  aerosol. 

The  global  total  surface  density  is  summarized  in  Table  2. 

Table  2.  Values  for  the  bin  sizes,  source  strengths  and  surface  area  for 
each  bin  size  for  A1203  particulates  from  SRM. 


Bin  Edges 

(/tm) 

Radius 

Dj  (jUm) 

Sj  (ton/yr) 

Global  Surface  Area,  A .  (10^  cm^) 

Baseline  Case  A  Case  B 

0.025-0.04 

0.03 

2.4 

3.50 

3.38 

1.09 

0.04-0.06 

0.05 

4.5 

4.22 

4.00 

1.94 

0.06-0.10 

0.08 

4.3 

2.64 

2.33 

1.60 

0.10-0.16 

0.13 

1.5 

0.62 

0.44 

0.39 

0.16-0.25 

0.20 

0.1 

0.035 

0.014 

0.014 

0.25-0.40 

0.32 

1.5 

0.22 

0.14 

0.14 

0.40-0.64 

0.51 

5.2 

0.48 

0.26 

0.26 

0.64-1.02 

0.81 

12.6 

0.73 

0.32 

0.32 

1.02-1.61 

1.28 

37.3 

1.33 

0.41 

0.41 

1.61-2.56 

2.03 

125.6 

2.84 

0.57 

0.57 

2.56-4.06 

3.23 

295.0 

4.28 

0.48 

0.48 

4.06-6.45 

5.12 

399.2 

3.74 

0.22 

0.22 

6.45-10.00 

8.13 

230.5 

1.43 

0.04 

0.04 

Total 

.... 

1120.0 

26.06 

12.60 

7.47 
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Section  4:  Input  from  Meteorites 


When  chunks  of  material  enter  the  mesosphere,  they  first  heated  and  ablated.  The 
exact  sequence  of  events  that  lead  to  the  production  of  micron  size  particles  is  rather 
complicated.  We  represent  this  incoming  material  by  a  flux  of  particles  deposited  at  the 
top  of  the  stratosphere.  Dr.  Lohn  (TRW)  estimates  that  at  the  top  of  the  stratosphere,  the 
mass  flux  of  one-micrometer  particle  from  meteorite  is  20-30  ton  per  year  with  an 
assumed  bulk  density  of  2.0  g/cm3,  this  corresponds  to  a  particle  flux  of  3.6  x  1018 
particles/year,  or  1.1  x  1011  particles/second.  We  took  the  upper  limit  of  the  data  as  input. 
The  input  particles  are  put  into  the  fifth  layer  from  top  of  the  model,  which  is 
approximately  from  55  km  to  60  km,  corresponding  to  the  top  of  the  stratosphere.  The 
source  is  distributed  evenly  around  the  globe  and  input  equally  everyday. 

For  the  Baseline  simulations  (no  sedimentation  and  no  removal  by  collision  with 
background  aerosols),  the  calculated  total  mass  is  1 13.25  ton.  Thus,  for  an  inert  tracer 
deposited  at  around  55-60  km,  the  residence  time  is  3.8  years.  The  corresponding  total 
surface  area  is  1.70  x  10  '  cm".  When  the  sedimentation  is  included  in  the  model,  the 

total  mass  is  31.5  ton.  The  calculated  residence  time  for  1pm  particles  with  pbulk  = 
2gm/cm  deposited  at  around  55-60  km  is  1.1  years.  The  total  surface  area  is  4.72  x  10 
cm".  We  did  not  show  the  results  from  the  case  with  collision  removal  because,  for 
particles  of  1  micrometer  size,  the  collision  removal  is  negligible  compared  with 
sedimentation  (see  Figure  6). 

The  latitude-altitude  distributions  for  4  seasons  of  the  zonally  mean  surface  area 
density  (10"  /An  1  cm  )  for  the  Baseline  simulation  and  Scenario  A  are  plotted  in  Figures 
1 1  and  12.  Because  of  the  smaller  input  (30  tons  versus  1120  tons)  the  total  surface  area 
density  in  the  meteor  source  case  is  more  than  one  order  of  magnitude  less  than  that  of 
SRM  case.  Compared  with  the  SRM  source  (Figures  12  and  10),  the  distribution  of  the 
total  surface  area  density  from  the  meteor  source  is  more  evenly  distributed  in  the  two 
hemispheres  because  the  uniform  source  distribution  assumption.  However,  there  are 
small  seasonal  variations  which  gives  higher  values  in  the  summer  hemisphere.  Similar 
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to  the  SRM  source  case,  the  lowest  values  of  the  total  surface  area  density  from  the 
meteor  source  occur  near  the  tropics.  The  latitude-longitude  distribution  of  the  surface 
area  density  (10'  jJm  /cm  )  at  20  km  from  Scenario  A  of  meteor  source  is  plotted  in 
Figure  13.  Compared  with  Figure  7,  Figure  13  shows  more  N-S  symmetry  and  also  the 
values  for  the  total  surface  area  density  are  much  smaller. 

Lohn  compared  the  debris  population/orbital  decay  methodology  presented  here  to 
the  results  of  Zolensky  [ 1989]  in  1994  (Lohn  et  al.,  [ 1994 ] ).  The  comparison  was 
reasonable;  we  calculated  particle  densities  in  the  stratosphere  (as  caused  by  deorbiting 
debris)  that  were  within  a  factor  of  three  of  the  measurements  of  Zolensky  [ 1989].  This 
work  was  described  previously  in  the  1994  report  (Lohn  et  al.,  [1994]). 
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Section  5:  Input  from  Orbital  Debris 


Dr.  Lohn  (TRW)  estimates  that  at  the  top  of  the  stratosphere,  the  mass  flux  of  one- 
micrometer  orbital  debris  objects  is  5-10  ton  per  year  with  the  density  of  4.7  g/cm  .  We 
took  the  upper  limit  of  the  data  as  input.  As  in  the  Meteor  case,  5.1  x  10  particles/year 
is  put  into  the  fifth  layer  from  top  of  the  model,  which  is  approximately  from  55  km  to  60 
km,  corresponding  to  the  top  of  the  stratosphere.  The  source  is  distributed  evenly  around 
the  globe  and  input  equally  everyday,  as  in  the  Meteor  case. 

For  the  Baseline  simulations  (no  sedimentation  and  no  removal  by  collision  with 
background  aerosols),  the  calculated  total  mass  of  the  orbital  debris  in  the  steady  state  is 
37.8  ton,  with  a  surface  area  of  2.41  x  101 1  cm2.  When  the  sedimentation  is  included  in 
the  model,  the  total  mass  of  the  orbital  debris  in  the  steady  state  is  6.1  ton.  The 
corresponding  residence  time  is  0.6  years.  Note  that  this  is  shorter  than  the  residence 

time  of  the  meteoritic  debris  because  pbulk  for  orbital  debris  is  assumed  to  be  a  factor  of 

2.3  larger.  The  total  surface  area  is  calculated  to  be  3.86  x  1010  cm2.  Again,  we  did  not 
include  collision  removal  in  the  background  sulfate  because  for  particles  of  1  micrometer 
size,  the  collision  removal  is  negligible  compared  with  sedimentation. 

The  latitude-altitude  distributions  for  4  seasons  of  the  zonally  mean  surface  area 
density  (10’  /Jm  1  cm  )  for  the  Baseline  case  and  Scenario  A  are  plotted  in  Figures  14  and 
15.  They  are  quite  similar  to  Figures  1 1  and  12,  because  the  particles  from  the  orbital 
debris  and  from  the  meteor  source  have  the  similar  source  distribution  and  same  sinks. 
The  only  difference  is  that  the  orbital  debris  has  a  smaller  source,  and  the  larger  bulk 
density  gives  even  smaller  number  density.  By  direct  scaling,  we  expect  the  surface  area 
to  be  a  factor  of  7  smaller  even  if  we  ignore  the  difference  in  residence  time. 
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Section  6:  Estimate  of  the  Effects  from  A1203  Particulate  on  Stratospheric  Ozone 
and  the  Sulfate  Layer 

A1203  particulate  may  impact  stratospheric  ozone  through  heterogeneous  reactions 
occurring  on  the  particle  surfaces.  The  reaction 

CIONO,  +  HCl  hno3  +  Cl2 

has  recently  been  measured  in  the  laboratory  by  Molina  et  al.,  [ 1997]  and  found  to 
proceed  at  a  reaction  probability  of  ~  0.2  on  aluminum  oxide.  A  similar  reaction  on  polar 
stratospheric  clouds  (PSCs)  and  sulfate  aerosol  is  believed  to  be  responsible  for  much  of 
the  Arctic  and  Antarctic  ozone  loss  (the  “ozone  hole”)  in  springtime.  Heterogeneous 
reaction  on  A1203  particles  would  not  be  limited  to  polar  regions  like  the  similar  reaction 
on  PSCs  and,  because  the  reaction  rate  is  not  temperature  dependent,  could  occur  at  all 
latitudes  and  seasons,  unlike  the  similar  sulfate  reaction.  The  above  reaction  converts 
C10N02  and  HCl,  both  moderately  long-lived  chlorine  reservoir  species,  into  Cl2,  a  very 
volatile  specie  that  dissociates  rapidly  in  the  presence  of  sunlight  into  atomic  Cl.  The 
chlorine  may  then  lead  to  ozone  removal  via  the  following  catalytic  cycle 

C10N02  +  HCl  -»  Cl2  +  HNO 3 
Cln  +  h  v  — >  2  Cl 
Cl  +  03  ->  CIO  +  02 
Cl  +  CH4  -»  HCl  +  CH3 
HNO:  +  hv  ->  OH  +  N02 
CIO  +  N02  — >  C10N02 

Net :  03  +  CH4  ->  02  +  OH  +  CH3 

The  impact  of  the  A1203  particulate  from  SRM  on  ozone  was  studied  with  a  2-D 
model  by  Jackman  et  al.,  [1998]  employing  the  Molina  et  al.,  [1997]  rate  for 
ClONCb+HCl  on  alumina  particles.  The  calculated  ozone  impact  is  extremely  small,  at 
most  0.06%  ozone  depletion  in  the  northern  polar  region  in  spring  and  a  global  average 
total  ozone  change  of  0.01%  depletion.  The  surface  area  density  of  A1203  used  by 
Jackman  et  al.,  [ 1998]  is  greater  than  our  calculated  surface  area  under  Case  B  by  factors 
of  2-10.  Therefore  we  can  state  that  ozone  depletion  on  the  global  scale  due  to  AHO3 
emissions  by  solid  rocket  motors  at  the  current  shuttle  launch  rate  is  less  than  0.04%  in 
the  northern  polar  region  and  ~  5  x  10"4  %  global  average. 
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Alumina  particles  from  SRM  may  affect  the  background  sulfate  aerosol  layer  in 
two  ways:  by  collision  of  A1203  particles  with  the  background  sulfate  layer,  increasing  its 
mass  and  surface  area,  and  by  acting  as  nuclei  for  H2SO4/H2O  condensation,  increasing 
the  number  of  sulfate  aerosol  particles  and  decreasing  their  mean  radius.  We  may  obtain 
an  upper  limit  estimate  of  the  first  effect  as  follows.  We  assume  that  those  A1203 
particles  that  collide  with  the  sulfate  particles  become  coated  with  sulfate  and  provide  the 
same  surface  area  for  sulfate  reactions.  Thus,  Case  A,  with  no  collision  removal,  and 
Case  B,  with  collision  removal  give  the  difference  in  surface  area.  The  maximum 
difference  is  4x1  O'4  pm2/cm3,  which  may  be  compared  with  background  sulfate  aerosol 
surface  areas  of  0. 1-1.0  pm2/ cm3.  Thus,  the  presence  of  A1203  particles  will  increase  the 
surface  area  of  sulfate  by  less  than  0.1%.  This  is  an  upper  limit  since  we  do  not 
compensate  for  decrease  in  sulfate  surface  area  for  the  sulfate  particles  that  collided  with 
the  A1203  particles. 

The  second  effect  can  be  evaluated  by  comparing  the  number  of  A1203  particles 
in  the  smallest  size  bin  (considered  to  be  the  only  size  likely  to  act  as  condensation 
nuclei)  with  the  total  number  of  sulfate  particles  per  unit  volume.  The  number  density  of 
A1203  particles  in  bin  1  is  plotted  in  Figure  16.  This  can  be  compared  with  sulfate 
particle  number  densities  of  1-20  particles  per  cubic  centimeter  for  those  particles  greater 
than  0.025  pm  radius.  We  conclude  that  the  impact  of  A1203  particles  as  condensation 
nuclei  for  sulfate  aerosol  is  probably  small  at  most  latitudes  and  altitudes,  but  has  the 
potential  to  increase  the  particle  number  density  and  surface  area  in  localized  regions 
where  number  density  is  otherwise  small.  It  should  be  noted  that  the  efficiency  of  A1203 
particles  to  act  as  condensation  nuclei  is  unknown,  so  assessing  this  effect  accurately  is 
not  possible  at  this  time. 

The  study  here  provides  a  good  indication  on  the  expected  effects  of  rocket 
launches.  A  more  realistic  evaluation  will  require  more  careful  treatments  of  the 
microphysical  processes  and  incorporation  of  more  realistic  launch  schedules. 
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Figure  1A:  The  simulated  steady  state  zonal  mean  mixing  ratio  of  Cl  (in  pptv)  due  to 
emission  from  9  space  shuttle  launches  annually.  It  is  assumed  that  the  first  launch 
occurs  on  day  40  of  the  year  with  one  launch  occurring  40  days  after.  Results  are  shown 
for  (a)  January,  (b)  April,  (c)  July,  and  (d)  October.  The  contour  levels  are  1,  3,  5,  7,  9, 
11  (pptv). 
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Figure  IB:  The  simulated  steady  state  zonal  mean  mixing  ratio  of  Cl  (in  pptv)  due  to 
emission  from  9  space  shuttle  launches  annually.  The  emissions  for  the  nine  launches  are 
put  into  the  model  at  a  uniform  rate  throughout  the  year.  Results  are  shown  for  (a) 
January,  (b)  April,  (c)  July,  and  (d)  October.  The  contour  levels  are  1,  3,  5,  7,  9  (pptv). 
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Figure  1C:  Model  calculated  steady  state  mixing  ratio  Cl  (in  pptv)  due  to  emission  from 
nine  shuttle  and  six  Titan  launches  per  year  as  repeated  in  Prather  et.al.  (1990b).  The 
total  Chlorine  input  is  725  tons/yr  compared  to  612  tons/yr  used  in  Figure  1A  and  IB. 
The  contour  levels  are  1,  3,  5, 1,9,  11  (pptv). 
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Figure  2:  The  sedimentation  velocity  (km/day)  of  the  AI2O3  particles  from  SRM  at  20°  N 
as  a  function  of  the  particle  size,  assuming  a  bulk  particle  density  of  1.7  gm/cc.  The 
dotted  curve  is  for  particles  at  40  km,  the  dashed  curve  for  30  km,  and  the  solid  curve  for 
20  km. 
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Figure  3:  The  collision  removal  rate  (yr1)  of  the  A1203  particles  from  SRM  at  20  km  as 
a  function  of  particle  size.  The  dotted  curve  is  for  particles  near  70°N,  the  dashed  curve 
is  for  particles  near  35°N,  and  the  solid  curve  is  for  particles  near  the  Equator. 
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Figure  4:  The  latitude-altitude  distribution  of  the  collision  removal  time  scale  (days)  of 
the  A1203  particles  in  the  first  bin  (radius=0. 025-0. 04  pm)  for  (a)  January,  (b)  April,  (c) 
July,  and  (d)  October. 
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Figure  5:  Same  as  Figure  4,  but  for  the  fourth  bin  (radius=0. 10-0. 16  urn). 
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Figure  6:  The  average  residence  time  (yr)  of  the  AI2O3  particles  from  SRM  as  a  function 
of  particle  size.  The  solid  curve  is  for  the  Baseline  Case,  the  dashed  curve  for  Case  A 
(sedimentation  only)  and  the  dotted  curve  for  Case  B  (sedimenation  and  collision 
removal).  Calculation  is  performed  using  the  size  bins  as  specified  in  Table  2. 
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Figure  7 :  The  latitude-longitude  distribution  of  the  total  surface  area  density 
(10  4  /Li  m2/cm3)  of  the  AI2O3  particles  from  SRM  at  20  km  in  January  for  Case  A 
(sedimentation  only).  The  contours  are  from  1  to  9  with  increment  of  1. 
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Figure  8:  The  latitude-altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10"  //  m7cm  )  of  the  A1203  particles  from  SRM  for  the  Baseline  Case  for  (a) 
January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours  are  from  1  to  9  with 
increment  of  1. 
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Figure  9:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10'  //  nT/cm  )  of  the  A1203  particles  from  SRM  for  Case  A  (sedimentation 
only)  for  (a)  January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours  are  from  1  to  9 
with  increment  of  1. 
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Figure  10:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10'  /./  m7cm  )  of  the  AI2O3  particles  from  SRM  for  Case  B  (sedimentation  and 
collision  removal)  for  (a)  January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours 
are  from  1  to  9  with  increment  of  1. 
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Figure  11:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10'  //  m  /cm  )  of  the  A1203  particles  from  meteor  source  for  the  Baseline  Case 
for  (a)  January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours  are  0.1  to  0.9  with 
increment  of  0.1. 
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Figure  12:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10  4  //  nr/cm3)  of  the  A1203  particles  from  meteor  for  Case  A  for  (a)  January, 
(b)  April,  (c)  July,  and  (d)  December.  The  contours  are  from  0.05  to  0.2  with  increment 
of  0.05. 
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Figure  13:  The  latitude-longitude  distribution  of  the  total  surface  area  density 

4  2  3 

(10  //  m"/cm  )  of  the  AI2O3  particles  from  meteor  source  at  20  km  in  January  for  Case  A 
(sedimentation  only).  The  contours  are  0.08,  0.1,  0.12,  0.15. 
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Figure  14:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  (10'  //  m  /cm  )  of  the  A1203  particles  from  orbit  debris  for  the  Baseline  Case  for 
(a)  January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours  are  from  0.02  to  0.1  with 
increment  of  0.02. 
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Figure  15:  The  latitude- altitude  distribution  of  the  zonally  averaged  total  surface  area 
density  ( 1 0  //  nT/cm  )  of  the  AI2O3  particles  from  orbit  debris  for  the  Case  A  for  (a) 
January,  (b)  April,  (c)  July,  and  (d)  December.  The  contours  are  0.005  and  0.01. 
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Figure  16:  The  latitude-altitude  distribution  of  the  zonally  averaged  number  density 
(10  2/cm3)  of  AI2O3  particles  in  bin  1  from  SRM  for  the  Baseline  case  for  (a)  January,  (b) 
April,  (c)  July,  and  (d)  October.  The  contours  are  from  0  to  2  in  increments  of  0.2. 
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